ORFEUS-II FAR-ULTRAVIOLET OBSERVATIONS OF 3C273: 
1. INTERSTELLAR AND INTERGALACTIC ABSORPTION LINESQ 
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ABSTRACT 

We present the first intermediate-resolution (A / 3000) spectrum of the bright quasi-stellar 
object 3C273 at wavelengths between 900 and 1200 A. Observations were performed with 

We detect 



the Berkeley spectrograph aboard the ORFEUS-II mission ( |Hurwitz et al. 1 



Lyman (3 counterparts to intergalactic Lyman a features identified by Morris et al. (1991) at 
cz = 19900, 1600, and 1000 km s _1 ; counterparts to other putative Lyman a clouds along the 
sight line are below our detection limit. The strengths of the two very low redshift Lyman (3 
features, which are believed to arise in Virgo intracluster gas, exceed preflight expectations 
( Wcymann et al. 1995| ), suggesting that the previous determination of the cloud parameters 
may underestimate the true column densities. A curve-of-growth analysis sets a minimum H I 
column density of 4 x 10 14 cm -2 for the 1600 km s _1 cloud. We find marginally significant 
evidence for Galactic H2 along the sight line, with a total column density of about 10 15 cm -2 . 
We detect the stronger interstellar O VI doublet member unambiguously; the weaker member is 
blended with other features. If the Doppler b value for O VI is comparable to that determined 



for N V (IScmbach, Savage, fc Tripp 1997| ) then the O VI column density is 7 ± 2 x 10 1 



significantly above the only previous estimate ( Davidscn 1993 ). The O VI/ N V ratio is about 
10, consistent with the low end of the range observed in the disk (compilation of Hurwitz 



Bowyer 1996). Additional interstellar species detected for the first time toward 3C273 (at 



modest statistical significance) include P II, Fe III, Ar I, and S III. 



Subject headings: ISM: atoms, galaxies: clusters: individual: Virgo, quasars: individual: 3C273 



1. OBSERVATIONS 

We present the first intermediate-resolution far-ultraviolet spectrum of the bright quasi-stellar object 



3C273 in the 900 - 1210 A band, emphasizing absorption features arising along the line of sight. Appenzeller 



st al. (199S ) discuss the intrinsic spectrum of 3C273. These data were collected with the Berkeley 



1 Based on the development and utilization of ORFEUS (Orbiting and Retrievable Far and Extreme Ultraviolet 
Spectrometers), a collaboration of the Institute for Astronomy and Astrophysics of the University of Tubingen, the Space 
Astrophysics Group of the University of California at Berkeley, and the Landessternwarte Heidelberg. 

2 Landessternwarte, University of Heidelberg 

institute for Astronomy and Astrophysics, University of Tubingen 
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spectrograph (|Hurwitz et al. 1998|) during the ORFEUS-SPAS II mission. The ORFEUS project and the 



ASTRO-SPAS platform are described in |Grewing et al. (199l| ) 



3C273 was observed five times, for a total of 10,797 s. All observations took place in late 1996, between 
330/21:41 and 338/16:33 (Day of Year/HH:MM, GMT). The target coordinates were a = 12 29 06.7, 6 = 
+02 03 09 (J2000). Absolute positioning of the 26" diameter ORFEUS entrance aperture is accurate to ± 



5". Extraction of the spectrum and subtraction of airglow follow the discussion in Hurwitz et al. (199S ). 
Statistical and systematic uncertainties associated with background and airglow subtraction are properly 
tracked. As a final step, we bin the data on 0.1 A centers. 

In Figure 1 we show the complete spectrum (solid line), the 1 a uncertainty associated with shot noise 
and detector flat-field effects (dotted line), and a continuum established by an automated fitting routine 
(dashed line). The automated fitting routine identifies absorption features deeper than a statistically 
determined threshold, replaces those data with a linear interpolation from nearby wavelengths, smoothes 
heavily, then iterates twice more at increasing sensitivity to absorption lines. The resulting continuum 
has the advantage of being determined objectively, and evidently tracks most of the true spectrum well at 
wavelengths greater than 1000 A. Like many such routines, however, this one yields unreliable results when 
strong features are closely spaced on a curved continuum (as in the region of absorption lines 3 and 4) and 
will tend to underestimate the true value when the spectrum contains a cluster of weak absorption features 
below the statistical threshold (note the weak dip in the fitted continuum around 1050 A). 

In Table 1 we list the interstellar absorption features and blends that are sufficiently strong and isolated 
for direct equivalent width measurement using the continuum shown in Figure 1 (unless otherwise noted). 
Statistical noise in the spectrum (lcr) corresponds to an unresolved absorption feature with equivalent 
width of about 0.040 A; continuum placement uncertainty corresponds to about half that value. The value 
in the fourth column includes the statistical uncertainty only. For species whose column densities have not 
previously been reported, we list the minimum column (if optically thin) based on the best fit equivalent 



width and the oscillator strcnths of Morton (1991) 



The identification of most of the strong interstellar features listed in Table 1 is secure ( Morton 1991 ). 
There are no obvious interstellar or intergalactic candidates for the unidentified line at 1008.5 A, which 
may be spurious. Although the observed features at 1066.8 and 1048.2 A presumably correspond to a pair 
of Ar I lines, their equivalent widths at first glance appear to be anomalous (the 1048.2 A line should be 
stronger). As noted, however, the longer wavelength feature may be partially blended with an H2 line 
at 1066.9 A. The shorter wavelength feature probably suffers from local continuum depression by nearby 
H2. Using the linear continuum shown in Figure 2, we estimate a best fit equivalent width closer to 0.13 
A. These systematic effects, in combination with the statistical errors, do much to reconcile the anomaly 
(especially if the features are on the flat part of the curve-of- growth). 

Below 1000 A, the increasing noise in the spectrum makes quantitative analysis of absorption features 
difficult, and the continuum fitting routine fails. Galactic features that are clearly detected include a blend 
from O I (and N III?) near 989 A, C III (977.02 A), Lyman S (972.54 A), and an N I blend near 953.8 A. 
Higher series Lyman lines are likely to be filled in by diffuse emission (we correct for diffuse H I emission 
only through Lyman S). An intriguing dip near 982 A may be associated with C III in the Virgo cluster, 
but this identification is speculative. Difficulties in continuum placement render the equivalent widths very 
uncertain; we attempt no quantitative analysis of absorption lines below 1000 A. 
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Table 1. Interstellar Lines and Equivalent Widths 



# W obs (A) W eq (A) a (A) ID Log (iV min ) (cm- 2 ) 



1 


1206 


.44 


0. 


.573 





.056 


Si III 1206.50 


13.4 


2 


1199 


.75 


0. 


.896 





,114 


N I blend 




3 a 


1193 


.56 


0. 


,510 





.050 


Si II 1193.29 




4 a 


1190 


.58 


0. 


,532 





.051 


Si II 1190.42 




5 


1152. 


.91 


0. 


.105 





.038 


P II 1152.82 


13.6 


6 


1144 


.92 


0. 


,340 





.054 


Fe II 1144.94 




7 


1143 


.44 


0, 


,137 





,040 


Fc II 1143.22 




8 


1134. 


.63 


1, 


,135 





.084 


N I blend 


15.1 


9 


1122 


.54 


0. 


,157 





.038 


Fe III 1122.53 


14.3 


10 


1121 


.99 


0. 


,149 





.038 


Fc II 1121.97 




11 


1096 


.94 


0, 


,202 





.047 


Fe II 1096.88 




12 


1084 


.19 


0. 


,515 





.051 


N II 1083.99, N II* blend 




13 b 


1066 


.77 


0. 


,135 





.047 


Ar I 1066.66 




14 


1063 


.20 


0. 


,327 





.049 


Fe II 1063.18 + H 2 




15 c 


1048 


.19 


0. 


,096 





,035 


Ar I 1048.22 


13.6 


16 


1039 


.18 


0. 


,433 





.044 


I 1039.23 




17 


1031 


.87 


0. 


,363 





.038 


VI 1031.93 




18 


1020. 


.75 


0. 


,129 





.040 


Si II 1020.70 




19 


1012. 


.53 


0. 


,083 





.036 


S III 1012.50 


14.4 


20 


1008. 


.51 


0. 


,165 





,039 


NOID 





a Continuum determined with 2nd order local polynomial fit. 
b Up to 0.02 A may be contributed by H 2 (J = 2). 
c Continuum may be depressed by adjacent weak lines. 
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INTERSTELLAR GAS 



The equivalent widths of features at wavelengths of overlap with the bandpass of the Hubble Space 



Telescope are consistent with previous measurements from that instrument (Morris et al. 1991, Savage 
1993). For several other interstellar species (such as Fe II), parameters inferred from the GHRS 



ct al. 



measurements ( Savage ct al. 1993] ) yield predictions for lines in the ORFEUS band that are consistent with 
the features observed. Our spectrum represents the first detection of P II, Fe III, Ar I, and S III along this 
sight line. The statistical significance of these features is modest, as can be seen in Table 1. The minimum 
columns inferred for these elements (and for N I and Si III, features of which have been observed previously 
but whose column densities have not been reported) are shown in Table 1 and are of moderate interest only. 
Our minimum column for Fe III is a factor of 2.5 below that established for Fe II ( Ravage ct al. 1993 ). The 
ratio of S III to S II ( Ravage et al. 1993| ) is at least 0.06 . 



Interesting interstellar species unique to the ORFEUS band include H2 and O VI, on which we now 
focus our attention. In Figure 2 we show a portion of the spectrum between 1046 and 1066 A. This region 
is comparatively clear of competing interstellar lines and provides a convenient hunting ground for features 
of H2. The solid line is the observed spectrum; the lowest dotted line is the 1 a error. Overlying the data is 
a synthetic spectrum consisting of a flat continuum plus interstellar features convolved with the instrument 
response (dashed line). Wherever possible we adopt column densities, velocities, and effective Doppler 



6- values of atomic species from Savage et al. (1993). An atomic-only model (omitted for clarity) does not 
reproduce the observed modulation near 1050 A nor the observed width of the the strong Fe II line at 
1063.2 A. 

The total H2 column density implied by our data is about 1 x 10 15 cm -2 , corresponding to a molecular 



fraction (/if 2 ) of about 2 x 10 . This value is not exceptional for a low-reddening sight line (Shull 



fc Beckwith 1982). The rotational excitation temperature T01 is poorly constrained but appears to be 
consistent with the range typical of interstellar clouds (50 - 80 K). 

In Figure 3 we show a portion of the spectrum between 1023 and 1041 A. The legend is as in Figure 2, 
except that we have not labeled a few weak features to avoid crowding. Features labeled Cld. 1 and Cld. 2 
arise in intergalactic gas. The discrepancy between their synthetic and observed spectra is discussed below. 
The strong galactic Lyman (3 line is well fit in the wings; at the core, the subtraction of bright diffuse 
emission introduces a significant uncertainty indicated by the large peak in the error tracing. The stronger 
O VI doublet member (1031.9 A) is resolved cleanly; the weaker member is blended but clearly present. To 
produce the fit shown in Figure 3, we found it necessary to set the column density of C II* near the upper 
limit from Savage et al. (1993). This modification to the C II* column does not significantly affect the 



previous determination that the cooling rate per nucleon is substantially below (about 1/6 of) the Galactic 



average (Savage et al. 1993) 



With a b value of about 35 km s 1 (estimated from the N V profile of Scmbach, Savage, & Tripp 1997), 
our best estimate of the O VI column density is 7 ± 2 x 10 14 cm~ 2 . Scmbach, Savage, & Tripp (1997) found 



an N V column of 8 x 10 13 cm~ 2 , yielding an O VI/ N V ratio of about 10. The ratio along disk sight lines 
varies from about 10 to 20 (compilation of Hurwitz fe Bowyer 1996| ). Ionization in this component of the 
hot galactic gas toward 3C273 therefore appears to be roughly similar to corresponding conditions in the 
disk. The fact that our measurement is near the low end of the disk range is qualitatively consistent with a 
scenario in which the gas is heated impulsively then cools as it moves away from the plane. Whether the 



comparatively high total column density is significantly nonrepresentative of the high latitude sky (Hurwitz 



fe Bowyer 1996) will be determined with certainty only after study of many additional sight lines. The 
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presence of Radio Loops I and IV (Berkhuijscn, Haslam, & Salter 1971) cautions at the very least that 
division of our projected column by the midplane density (Jenkins 1978) should not be taken as a reliable 
estimate of the typical scale height. However, the fact that this first clear measurement of columns for both 
O VI and N V yields a value comparable to that in the disk bodes well for future far-ultraviolet missions 
such as FUSE. The previous measurement of O VI toward 3C273 ( Davidson 1993| ) had set only a lower limit 
on the column density, corresponding to an O VI/ N V ratio substantially below the range of disk values. 



3. INTERGALACTIC GAS 



Sixteen otherwise unidentified absorption features in the GHRS spectrum were attributed to 



intergalactic Lyman a clouds by Morris et al. (1991). Most of these are sufficiently weak that their higher 
Lyman series lines would not be detectable in our spectrum. We do find probable Lyman (3 counterparts 



to the strongest proposed Lyman a forest features, however. Morris et al. (1991 ) propose a cloud at cz 
= 19900 km s , from its column density and b value we predict a Lyman f3 equivalent width of 0.080 
A at 1093.93 A. We find an otherwise unidentified feature within 0.2 A of the expected position with an 
equivalent width of 0.060 ± 0.038 A confirming, at least to 1.6 a, the intergalactic H I origin. 

Of greater interest are a strong pair of features attributed to comparatively low redshift H I gas in the 
Virgo supercluster. Their Lyman a lines appear near +1000 and +1600 km s _1 (Weymann et al. 1995). 
Lyman (3 counterparts are expected at 1029.21 and 1031.16 A. We detect absorption features at 1029.11 
and 1031.14 A; however, the equivalent widths are stronger than expected. Weymann et al. (1995[ ) employ 
profile fitting to constrain column densities and Doppler b values. Based on those parameters, the Lyman 
(3 equivalent widths are predicted to be 0.098 and 0.10 A, respectively. We observe 0.145 (± 0.037) and 
0.241 (± 0.032) A. The disagreement is seen clearly in Figure 3, where the synthetic spectrum relies on the 
Weymann ct al. (1995| ) parameters. 



The discrepancy is particularly pronounced for the +1600 km s 1 cloud. Even with the parameters 
proposed by Weymann et al. the Lyman a line is saturated (tq = 3.7), so determination of a column 



density based solely on profile fitting of that feature may be affected by unresolved velocity structure 
and/or uncertainties in the instrumental profile. We have analyzed the GHRS spectrum to extract the 
equivalent width of the Lyman a line. In combination with the Lyman (3 line from the ORFEUS spectrum, 
we establish constraints on the logarithmic column density and effective b value for the +1600 km s _1 cloud 



using a curve-of-growth technique. We show the results, and the parameters proposed by Weymann ct al 



in Figure 4. Our results suggest that if the feature arises in a single cloud, its column density is higher by at 
least a factor of four compared to the value suggested by Weymann et al, . The curve-of-growth results are 



consistent with b values as low as 12 km s _1 and column densities in excess of 10 18 cm -2 . Although this 
extremum of parameter space is probably inconsistent with the Lyman a profile fitting, it is not difficult to 
construct a multiple-cloud system that closely approximates the theoretical Voigt profile inferred from the 



Weymann et al. parameters while hiding a great deal of gas in the saturated core. 



1 



van Gorkom et al. (1993) searched for 21 cm emission from Virgo gas near the velocities of the 
absorbing clouds, detecting no H I to a limit of about 2.8 x 10 19 cm -2 . By adopting a relationship between 
the volume density of clouds with columns near the radio limit and of clouds with columns near those of 



the absorption line systems (Tytler 1987), the radio survey constrains the minimum size of the absorbing 



systems. The minimum permitted cloud size grows slowly with column density (r c i ou d oc N ), so our 



results increase the minimum cloud size from about 4 kpc (van Gorkom et al. 1993) only to about 6 kpc, at 
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least for the smallest columns within our 95% confidence interval. 

Other interstellar species that might conceivably be blended with the +1600 km s _1 cloud include CI I 
at 1030.88 A, or H 2 (</ = 3) at 1030.89 A. Based on the nondetection of other features from these species, 
we conservatively limit the equivalent width of their combined contribution to no more than 0.050 A. A 
contribution from Galactic O VI at high negative velocity is more difficult to rule out formally from our 
spectrum. Neither N V nor C IV appears at high negative velocity, however ( Sembach, Savage, fc Tripp| 
1997b. 



CONCLUSIONS 



We observed the bright quasi-stellar object 3C273 for 10,797 s with the Berkeley spectrograph aboard 
the ORFEUS-II mission ( Hurwitz et al. 1998 ). The resulting spectrum offers intermediate spectral 
resolution (A / 3000 FWHM) and a la noise corresponding to an absorption feature with equivalent width 
of about 0.040 A. 

The spectrum reveals, at modest statistical significance, a variety of previously undetected low 
ionization interstellar species, including P II, Fe III, Ar I, and S III. There is marginally significant evidence 
for H2, with a total column density of about 10 15 cm~ 2 . The J = and J = 1 rotational levels seem 
comparably populated, consistent with an excitation temperature Xbi typical for interstellar clouds (50 - 
80 K). 

We detect the stronger interstellar O VI doublet member unambiguously; the weaker member is 



blended. If the Doppler b value for O VI is comparable to that determined for N V (Sembach, Savage, & 
Tripp 1997|), the O VI column density is 7 ± 2 x 10 14 cur 2 . The O VI/ N V ratio is about 10, consistent 



with the low end of the range observed in the disk (compilation of Hurwitz & Bowyer 1996). If impulsively 
heated gas cools as it moves away from the Galactic disk, the O VI/ N V ratio might be expected to decrease 
(compared to the disk value) along halo sight lines. 



We detect Lyman (3 counterparts to intergalactic Lyman a features identified by Morris et al. (1991) 
at cz = 19900, 1600, and 1000 km s _1 . Counterparts to other putative Lyman a clouds along the sight 
line are below our detection limit. The strengths of the two very low redshift Lyman (3 features, which 
are believed to arise in Virgo intracluster gas, exceed preflight expectations ( Weymann et al. 1995 ). A 
curve-of-growth analysis sets a minimum H I column density of 4 x 10 14 cm -2 for the 1600 km s _1 cloud. 



Following van Gorkom et al. (1993), our revised column density and nondetection of 21 cm radio emission 
discussed in that work can set a lower limit of about 6 kpc to the size of the Virgo cluster clouds. 



We acknowledge our colleagues on the ORFEUS team and the many NASA and DARA personnel who 
helped make the ORFEUS-II mission successful. This work is supported by NASA grant NAG5-696. 
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Fig. 1. — Observed spectrum of 3C 273 (solid line), 1 a uncertainty (dotted line), continuum fit longward of 
1000 A(dashed line). Numbered interstellar absorption features are identified in Table 1. The modulation 
near 1069 A is an artifact associated with second-order He I diffuse emission. "IG" indicates intergalactic 
Lyman [5 features. 
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Fig. 2. — Observed spectrum of 3C 273 (solid line) with 1 a uncertainty (lower dotted line). Synthetic 
spectra overlying the data are based on a flat continuum plus interstellar atomic and H 2 features (dashed 
line). H 2 features are labeled with J level of lower state. 
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Fig. 3. — Observed spectrum of 3C 273 (solid line) with 1 a uncertainty (dotted line). Synthetic spectrum 
overlying the data is based on a flat continuum plus interstellar atomic and H 2 features and intergalactic 
features. H2 features are labeled with J level of lower state. Features labeled "Cld. 1" and "Cld. 2" are 
Lyman (3 clouds in the Virgo cluster. 
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Fig. 4. — Constraints on logarithmic column density and effective Dopplcr b value for the 1600 km s _1 
H I cloud. The cross centered at (14.22, 34) is based on Lyman a profile fitting (Weymann et al. 1995). 
The best-fit, 68% and 95% confidence intervals based on the equivalent widths of Lyman a (GHRS) and f3 
(ORFEUS) are also shown. 



